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Preface 



This book deals primarily with fundamental aspects of engineering vibrations 
within the framework of the linear theory. Although it is true that in practical 
cases it is sometimes not easy to distinguish between linear and nonlinear 
phenomena, the basic assumption throughout this text is that the principle 
of superposition holds. 

Without claim of completeness, the authors’ intention has been to discuss 
a number of important topics of the subject matter by bringing together, in 
book form, a central set of ideas, concepts and methods which form the 
common background of real-world applications in disciplines such as 
structural dynamics, mechanical, aerospace, automotive and civil engineering, 
to name a few. 

In all, the authors claim no originality for the material presented. However, 
we feel that a book such as this one can be published at the end of the 1990s 
because, while it is true that the general theory of linear vibrations is well 
established (Lord Rayleigh’s book Theory of Sound is about a century old), 
this by no means implies that the subject is ‘closed’ and outside the mainstream 
of ongoing research. In fact, on the one hand, the general approach to the 
subject has significantly changed in the last 30 years or so. On the other 
hand, the increasing complexity of practical problems puts ever higher 
demands on the professional vibration engineer who, in turn, should acquire 
a good knowledge in a number of disciplines which are often perceived as 
distinct and separate fields. 

Also, in this regard, it should be considered that the computer revolution 
of recent years, together with the development of sophisticated algorithms 
and fully automated testing systems, provide the analyst with computation 
capabilities that were unimaginable only a few decades ago. This state of 
affairs, however — despite the obvious advantages — may simply lead to 
confusion and/or erroneous results if the phenomena under study and the 
basic assumptions of the analysis procedures are not clearly understood. 

The book is divided into two parts. Part I (Chapters 1 to 12) has been 
written by Paolo L.Gatti and is concerned with the theory and methods of 
linear engineering vibrations, presenting the topics in order of increasing 
difficultly — from single-degree-of-freedom systems to random vibrations and 
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stochastic processes — and also including a number of worked examples in 
every chapter. Within this part, the first three chapters consider, respectively, 
some basic definitions and concepts to be used throughout the book (Chapter 
1), a number of important aspects of mathematical nature (Chapter 2) and 
a concise treatment of analytical mechanics (Chapter 3). In a first reading, 
if the reader is already at ease with Fourier series, Fourier and Laplace 
transforms, Chapter 2 can be skipped without loss of continuity. However, 
it is assumed that the reader is familiar with fundamental university calculus, 
matrix analysis (although Appendix A is dedicated to this topic) and with 
some basic notions of probability and statistics. 

Part II (Chapters 13 to 15) has been written by Vittorio Ferrari and deals 
with the measurement of vibrations by means of modern electronic 
instrumentation. The reason why this practical aspect of the subject has been 
included as a complement to Part I lies in the importance — which is sometimes 
overlooked — of performing valid measurements as a fundamental requirement 
for any further analysis. Ultimately, any method of analysis, no matter how 
sophisticated, is limited by the quality of the raw measurement data at its 
input, and there is no way to fix a set of poor measurements. The quality of 
measurement data, in turn, depends to a large extent on how properly the 
available instrumentation is used to set up a measuring chain in which each 
significant source of error is recognized and minimized. This is especially 
important in the professional world where, due to a number of reasons such 
as limited budgets, strict deadlines in the presentation of results and/or real 
operating difficulties, the experimenter is seldom given a second chance. 

The choice of the topics covered in Part II and the approach used in the 
exposition reflect the author’s intention of focusing the attention on basic 
concepts and principles, rather than presenting a set of notions or getting 
too much involved in inessential technological details. The aim and hope is, 
first, to help the reader — who is only assumed to have a knowledge of basic 
electronics — in developing an understanding of the essential aspects related 
to the measurement of vibrations, from the proper choice of transducers and 
instruments to their correct use, and, second, to provide the experimenter 
with guidelines and advice on how to accomplish the measurement task. 

Finally, it is possible that this book, despite the attention paid to reviewing 
all the material, will contain errors, omissions, oversights and/or misprints. 
We will be grateful to readers who spot any of the above or who have any 
comment for improving the book. Any suggestion will be received and 
considered. 

Milan 1998 

Paolo Luciano Gatti, 
Vittorio Ferrari 

Email addresses: 
pljgatti@tin.it 
ferrari@bsing.ing.unibs.it 
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Paolo L.Gatti 



Copyright © 2003 Taylor & Francis Group LLC 




1 



Review of some fundamentals 



1.1 Introduction 

It is now known from basic physics that force and motion are strictly 
connected and are, by nature, inseparable. This is not an obvious fact; it 
has taken almost two millennia of civilized human history and the effort 
of many great minds to understand. At present, it is the starting point of 
almost every branch of known physics and engineering. One of these 
branches is dynamics: the study that relates the motion of physical bodies 
to the forces acting on them. Within certain limitations, this is the realm of 
Newton’s laws, in the framework of the theory that is generally referred to 
as classical physics. 

Mathematically, the fact that force causes a change in the motion of a 
body is written 

d 

F = — (mv) (1.1) 

dt 

This is Newton’s second law which defines the unit of force once the 
fundamental units of mass and distance are given. 

An important part of dynamics is the analysis and prediction of vibratory 
motion of physical systems, in which the system under study oscillates about 
a stable equilibrium position as a consequence of a perturbing disturbance 
which, in turn, starts the motion by displacing the system from such a position. 

This type of behaviour and many of its aspects — wanted or unwanted — 
is common everyday experience for all of us and is the subject of this book. 
However, it must be clear from the outset that we will only restrict our 
attention to ‘linear vibrations’ or, more precisely, to situations in which 
vibrating systems can be modelled as ‘linear’ so that the principle of 
superposition applies. Future sections of this chapter and future chapters 
will clarify this point in stricter detail. 
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1.2 The role of modelling (linear and nonlinear, discrete 
and continuous systems, deterministic and random 
data) 

In order to achieve useful results, one must resort to models. This is true in 
general and applies also to all the cases of our concern. Whether these models 
be mathematical or nonmathematical in nature, they always represent an 
idealization of the actual physical system, since they are based on a set of 
assumptions and have limits of validity that must be specified at some point 
of the investigation. So, for the same system it is possible to construct a 
number of models, the ‘best’ being the simplest one that retains all the essential 
features of the actual system under study. 

Generally speaking, the modelling process can be viewed as the first step 
involved in the analysis of problems in science and engineering: the so-called 
‘posing of the problem’. Many times this first step presents considerable 
difficulties and plays a key role to the success or failure of all subsequent 
procedures of symbolic calculations and statement of the answer. With this 
in mind, we can classify oscillatory systems according to a few basic criteria. 
They are not absolute but turn out to be useful in different situations and 
for different types of vibrations. 

First, according to their behaviour, systems can be linear or nonlinear. 
Formally, linear systems obey differential equations where the dependent 
variables appear to the first power only, and without their cross products; 
the system is nonlinear if there are powers greater than one, or fractional 
powers. When the equation contains terms in which the independent variable 
appears to powers higher than one or to fractional powers, the equation 
(and thus the physical system that the equation describes) is with variable 
coefficients and not necessarily nonlinear. The fundamental fact is that for 
linear system the principle of superposition applies: the response to different 
excitations can be added linearly and homogeneously. In equation form, if 
f(x) is the output to an input x, then the system is linear if for any two inputs 
x, and x 2 , and any constant a, 



f{xi +X 2 ) =f(x\) + fixz) 


(additive property) 


(1.2) 


f(ax) = af{x) 


(homogeneous property) 


(1.3) 



The distinction is not an intrinsic property of the system but depends on 
the range of operation: for large amplitudes of vibration geometrical 
nonlinearity ensues and — in structural dynamics — when the stress-strain 
relationship is not linear material nonlinearity must be taken into account. 

Our attention will be focused on linear systems. For nonlinear ones it is 
the author’s belief that there is no comprehensive theory (it may be argued 
that this could be their attraction), and the interested reader should refer to 
specific literature. 
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Second, according to the physical characteristics — called parameters — 
systems can be continuous or discrete. Real systems are generally continuous 
since their mass and elasticity are distributed. In many cases, however, it is 
useful and advisable to replace the distributed characteristics with discrete 
ones; this simplifies the analysis because ordinary differential equations for 
discrete systems are easier to solve than the partial differential equations 
that describe continuous ones. 

Discrete-parameter systems have a finite number of degrees of freedom, 
i.e. only a finite number of independent coordinates is necessary to define 
their motion. The well-known finite-element method, for example, is in 
essence a discretization procedure that retains aspects of either continuous 
and discrete systems and exploits the calculation capabilities of high-speed 
digital computers. Whatever discretization method is used, one advantage is 
the possibility to improve the accuracy of the analysis by increasing the 
number of degrees of freedom. 

Also in this case, the distinction is more apparent than real; continuous systems 
can be seen as limiting cases of discrete ones and the connection of one 
formulation to the other is very close. However, a detailed treatment of continuous 
systems probably gives more physical insight in understanding the ‘standing- 
wave-travelling-wave duality’, intrinsic in every vibration phenomenon. 

Third, in studying the response of a system to a given excitation, sometimes 
the type of excitation dictates the analysis procedure rather than the system 
itself. From this point of view, a classification between deterministic and 
random (or stochastic or nondeterministic) data can be made. Broadly 
speaking, deterministic data are those that can be described by an explicit 
mathematical relationship, while there is no way to predict an exact value at 
a future instant of time for random data. In practice, the ability to reproduce 
the data by a controlled experiment is the general criterion to distinguish 
between the two. With random data each observation will be unique and 
their description is made only in terms of statistical statements. 

1.3 Some definitions and methods 

As stated in the introduction, the particular behaviour of a particle, a body or 
a complex system that moves about an equilibrium position is called oscillatory 
motion. It is natural to try to describe such a particle, body or system using an 
appropriate function of time x(t). The physical meaning of x(t) depends on 
the scope of the investigation and, as often happens in practice, on the available 
measuring instrumentation: it might be displacement, velocity, acceleration, 
stress or strain in structural dynamics, pressure or density in acoustics, current 
or voltage in electronics or any other quantity that varies with time. 

A function that repeats itself exactly after certain intervals of time is called 
periodic. The simplest case of periodic motion is the harmonic (or sinusoidal) 
that can be defined mathematically by a sine or cosine function: 
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Jc(i) = X cos (ut - 0) 



(1.4) 



where: 

X is the maximum, or peak amplitude (in the appropriate units); 

(cot-0) is the phase angle (in radians); 
ft) is the angular frequency (in rad/s); 

9 is the initial phase angle (in radians), which depends on the choice of 
the time origin and can be taken equal to zero if there is no relative 
reference to other sinusoidal functions. 



The time between two identical conditions of motion is the period T. It is 
measured in seconds and is the inverse of the frequency v whose unit is the 
hertz (Hz, with dimensions of s -1 ) and, in turn, represents the number of 
cycles per unit time. The following relations hold: 



w — 2 TCV 



(1.5) 



1 2tt 
v to 



( 1 . 6 ) 



A plot of eq (1.4), amplitude versus time, is shown in Fig. 1.1 where the 
peak amplitude is equal to unity. 

A vector x of modulus X that rotates with angular velocity coin the xy plane 
is a useful representation of sinusoidal motion: x(t) is now the instantaneous 
projection of x on the x-axis (on the y-axis for a sine function). 




Fig. 1.1 Cosine harmonic oscillation of unit amplitude. 
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Other representations are possible, each one with its own particular 
advantages; however, the use of complex numbers for oscillatory quantities 
gives probably the most elegant and compact way of dealing with the 
problem. 

Recalling from basic calculus the Euler equations 

e~‘ z = cos z — i sin z 

1 . 

cos 2 = - (e ,z + e w ) 

2 

1 

sinz = — {e* z — e ,z ) 

2 i (1.7) 

where i is the imaginary unit and (i = 'f— 1) e=2. 71828. ..is the well-known 
basis of Naperian logarithms, an oscillatory quantity can be conveniently 
written as the complex number 

x(t) - Ce~^ (1-8) 

where C is the complex amplitude, i.e. a complex number that contains 
both magnitude and phase information and can be written as (a+ib) or Xe w 
j C | = X = y/a 2 + b 1 with magnitude and phase angle 9, where and 9=b/a, 

cos 9=a/X and sin 9=b/X. 

The number C* = a — ib is the complex conjugate of C and the square of 
the magnitude is given by | C | 2 = CC* = X 2 =a 2 + b z . 

The idea of eq (1.8) — called the phasor representation — is the temporary 
replacement of a real physical quantity by a complex number for purposes 
of calculation; the usual convention is to assign physical significance only to 
the real part of eq (1.8), so that the oscillatory quantity x(t) can be expressed 
in any of the four ways 

x(t) = a cos( U3t ) + b sin( wt) = X cos( ut — 9) ( 1 • 9 ) 

or 

xtf) = Ce“^=Xe- ) ' M - ()| (1-10) 

where only the real part is taken of the expressions in eq (1.10). 

A little attention must be paid when we deal with the energy associated 
with these oscillatory motions. The various forms of energy (energy, energy 
density, power or intensity) depend quadratically on the vibration amplitudes, 
and since Re(x 2 ) ^ [Re(x)] 2 we need to take the real part first and then square 
to find the energy. 

Furthermore, it is often useful to know the time-averaged energy or 
power and there is a convenient way to extract this value in the general 
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case. Suppose we have the two physical quantities expressed in the real 
form of eq (1.4) 

x\ (() = Xj cosher - 81 ) 

*2(0 = X2 COS(U2t - 82) 

It is easy to show that {* 1 X 2 }) he. the average value of the product * 1 X 2 , is 
different from zero only if fjJ\ = UJi and in this case we get 

1 f 7 ' 1 

(x\xi) = — xi(t)xz(t)dt = - XiX 2 cos( 8 i - 82 ) (1.11) 

T Jo 2 

where T = lirjw and the factor 1/2 comes from the result (cos 2 ) = 1/2. 

If we want to use phasors and represent the physical quantities as 

* ! (t) = Xie- )W -* 11 
x 2 (t} = X 2 e- il ^ l - ei] 

we see that in order to get the correct result of eq (1.11) we must calculate 
the quantity 

^Refoxf) (1.12) 

In the particular case of X\ —Xi^X (where these terms are expressed in 
the form of eq (1.8)), our convention says that the average value of x squared 
is given by 

jRefxx*) = \ CC* = \\ C\ z =jX 2 =\{a 2 + b 2 ) (1.13) 

Phasors are very useful for representing oscillating quantities obeying linear 
equations; other authors (especially in electrical engineering books) use the 
letter / instead of i and & m instead of g - and some other authors use the 

positive exponential notation e im . Since we mean to take the real part of the 
result, the choice is but a convention; any expression is fine as long as we 
are consistent. The negative exponential is perhaps more satisfactory when 
dealing with wave motion, but in any case it is possible to change the formulas 
to the electrical engineering notation by replacing every i with —j. 

Periodic functions in general are defined by the relation 

x(t) = x(f + nT) n = 1,2, 3, ... (1.14) 

and will be considered in subsequent chapters where the powerful tool of 
Fourier analysis will be introduced. 
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1.3.1 The phenomenon of beats 

Let us consider what happens when we add two sinusoidal functions of slightly 
different frequencies w l and iv 2 , with ^2 — Wi + e and e being a small quantity 
compared to w 1 and u > 2 . In phasor notation, assuming for simplicity equal 
magnitude and zero initial phase for both oscillations and x 2 , we get 

{t)+x 2 {t) = Xe-'™' ‘ + Xe~ k “' lt ( 1 ■ ■ 1 ■ 5 ) 



that can be written as 



/(us - ui\ )tj2 _[_ e ~ Huii - ui\)t]2 j 



with a real part given by 

2X cos(wf)cos(w fltt f) (1.16) 

where w ai , = (wt +wi)/2 and S = (tC 2 — We can see eq (1.16) as an 

oscillation of frequency w av and a time-dependent amplitude 2X cos (wt). A 
graph of this quantity is shown in Figs 1.2 and 1.3 where wj — 8,0 rad/s and 
u ?2 — wi = 0.6 (Fig. 1.2) and 1.0 (Fig. 1.3), respectively. 

Physically, the two original waves remain nearly in phase for a certain 
time and reinforce each other; after a while, however, the crests of the first 
wave correspond to the troughs of the other and they practically cancel out. 
This pattern repeats on and on and the result is the phenomenon of beats 
illustrated in Figs 1.2 and 1.3. Maximum amplitude occurs when fit = mv 
(n= 0, 1, 2,...), that is, every tt/ 2’ seconds. Therefore, the frequency of the 




Fig. 1.2 Beat phenomenon (u >2 — u>] = 0.6), 
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Fig. 1.3 Beat phenomenon (wj — u!\ = 1.0). 



beats is (w/tt) = vi — i'i, equal to the difference in frequency between the two 
original signals. 

When the signals have unequal amplitudes (say A and B) the total 
amplitude does not become zero, it oscillates between A+B and A-B, but the 
general pattern can still be easily identified. 



CtitifMna grami* (h* 32 dlr.x} 




MW 14.M Ifr.SW 1*19$ 3*1*4 



Fig. 1.4 Beat phenomenon measured on an Italian belltower (oscillation of largest 
bell). 
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In acoustics, for example, beats are heard as a slow rise and fall in sound 
intensity (at the beat frequency l 'l ~ v 2) when two notes are slightly out of 
tune. Many musicians exploit this phenomenon for tuning purposes, they 
play the two notes simultaneously and tune one until the beats disappear. 

Figure 1.4 illustrates a totally different situation. It is an actual vibration 
measurement performed on an ancient belltower in Northern Italy during 
the oscillation of the biggest bell. The measurement was made at about two- 
thirds of the total height of about 50 m on the body of the tower in the 
transverse direction. There is a clear beat phenomenon between the force 
imposed on the structure by the oscillating bell (at about 0.8 Hz, with little 
variations of a few percent) and the first flexural mode of the tower (0.83 
Hz). Several measurements were made and the beat frequency was shown to 
vary between 0.03 and 0.07 Hz, indicating a situation close to resonance. 
This latter concept, together with the concepts of forced oscillations and 
modes of a vibrating system will be considered in future chapters. 

1.3.2 Displacement, velocity and acceleration 

If the oscillating quantity x(t) in eq (1.8) is a displacement and we recall the 
usual definitions of velocity and acceleration 

dx(t) 

v(t)=x(t) = —— 

dt 

a(t)=x{t) = 

we get from the phasor representation 

w(f) = -iuCe-** = wGr i(wt+ ^ 2) 

a<t) = ^ 2 Ce-^‘ + (1 ' 17) 

since the phase angle of -i is n/2 and the phase angle of -1 is n. The velocity 
leads the displacement of 90°, the acceleration leads the velocity of 90° and 
all three rotate clockwise in the Arland-Gauss plane (abscissa=real part, 
ordinate=imaginary part) as time passes. Moreover, from eqs (1.17) we note 
that the maximum velocity amplitude is Vim* — wC, while the maximum 
acceleration amplitude is = w 2 C. 

In theory it should not really matter which one of these three quantities 
is measured; the necessary information and frequency content of a signal is 
the same whether displacement, velocity or acceleration is considered and 
any one quantity can be obtained from any other one by integration or 
differentiation. However, physical considerations on the nature of vibrations 




Copyright © 2003 Taylor & Francis Group LLC 



themselves and on the electronic sensors and transducers used to measure 
them somehow make one parameter preferred over the others. 

Physically, it will be seen that displacement measurements give most weight 
to low-frequency components and, conversely, acceleration give most weight 
to high-frequency components. So, the frequency range of the expected signals 
is a first aspect to consider; when a wide-band signal is expected, velocity is 
the appropriate parameter to select because it weights equally low- and high- 
frequency components. Furthermore, velocity (rms values, see next section), 
being directly related to the kinetic energy, is preferred to quantify the severity, 
i.e. the destructive effect, of vibration. 

On the other hand, acceleration sensitive transducers (accelerometers) are 
commonly used in practice because of their versatility: small physical 
dimensions, wide frequency and dynamic ranges, easy commercial availability 
and the fact that analogue electronic integration is more reliable than 
electronic differentiation are important characteristics that very often make 
acceleration the measured parameter. All these aspects play an important 
role but, primarily, it must be the final scope and aim of the investigation 
that dictates the choice to make for the particular problem at hand. 

Let us make some heuristic considerations from the practical point of 
view of the measurement engineer. Suppose we have to measure a vibration 
phenomenon which occurs at about v=l Hz with an expected displacement 
amplitude (eq (1.17)) of C=±l mm. It is not difficult to find on the market 
a cheap displacement sensor with, say, a total range of 10 mm and a sensitivity 
of 0.5 V/mm. In our situation, such a sensor would produce an output signal 
of 1 V, meaning a good signal-to-noise ratio in most practical situations. On 
the other hand, the acceleration amplitude in the above conditions is about 
(2tt0 2 C = (39.48)(2 x 10“ 3 ) = 7.9 x 10~ 2 m/s 2 ^8.1 x 10“ 3 g so that a 
standard general-purpose accelerometer with a sensitivity of, say, 100 mV/g 
would produce an output signal of 0.81 mV — 8,1 x 10 -4 V, which is much 
less favourable from a signal-to-noise ratio viewpoint. 

By contrast — for example in heavy machinery — forces applied to massive 
elements generally result in small displacements which occur at relatively 
high frequencies. So, for purpose of illustration, suppose that a machinery 
element vibrates at about 100 Hz with a displacement amplitude of ±0.05 
mm. The easiest solution in this case would be an acceleration measurement 
because the acceleration amplitude is now 39.5 m/s 2 = 4g, so that a general 
purpose (and relatively unexpensive) 100 mV/g accelerometer would produce 
an excellent peak-to-peak signal of about 400 mV. In order to measure such 
small displacements at those values of frequency, we would probably have 
to resort to more expensive optical sensors. 

1.3.3 Quantification of vibration level and the decibel scale 

The most useful descriptive quantity — which is related to the power content 
of the of the vibration and takes the time history into account — is the root 
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mean square value (x rm J, defined by 



^rmS — 



\ 



lim — 

T^-k, J J 



x 2 (t)dt 



(1.18) 



For sinusoidal motion it is easily seen that x rms = Xj-Jl = 0.707X. In the 
general case X/x rms is called the the form factor and gives some indication of 
the waveshape under study when impulsive components are present or the 
waveform is extremely jagged. It is left as an easy exercise to show that for 
a triangular wave (see Fig. 1.5) x rms =0.577X. 

In common sound and vibration analysis, amplitudes may vary over wide 
ranges (the so-called dynamic range) that span more than two or three 
decades. Since the dynamic range of electronic instrumentation is limited 
and the graphical presentation of such signals can be impractical on a linear 
scale, the logarithmic decibel (dB) scale is widely used. 

By definition, two quantities differ by one bel if one is 10 times ( 10 1 ) 
greater than the other, two bels if one is 100 times ( 10 2 ) greater than the 
other and so on. One tenth of a bel is the decibel and the dB level (L) of a 
quantity x, with respect to a reference value x 0 , is given by 



L(dB) = 10 log 10 




(1.19) 
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So, for example, 3 dB represent a doubling of the relative quantity, i.e. 
x = 2.0xo.Decibels, like radians, are dimensionless; they are not ‘units’ in the 
usual sense and consistency dictates that the reference value must be 
universally accepted. When this is not the case, the reference value must be 
clearly specified. Sound intensity levels (SIL), in acoustics, are given by eq 
(1.19) and the reference intensity value is /o = 10“ 12 W/m 2 , 

A more commonly used version of this logarithmic scale requires the 
proportionality factor 20 instead of 10, and is related to the fact — already 
mentioned in the preceding paragraph — that energies, powers, rms values, 
etc. depend quadratically on the amplitudes. 

We have 

L(dB) = 20 log 10 




In acoustics, eq (1.20) defines the sound pressure level (SPL) where the reference 
value (in air) is the pressure po = 2 x 10 -J Pa. In vibration measurements one 
can have a displacement level L d , a velocity level L,„ an acceleration level L a or 
a force level L f . These are obtained, respectively, by the relations 

L d = 20\og lo {d/d 0 ) 

L v = 201og 10 (tVt> 0 ) 

L a = 20 !og 10 (<i/tf 0 ) 

Lf — 20 log 10 (/" /fa) 

The reference values are not universally accepted, but usually we have 

do = 10 -11 m 
v 0 = 10 -9 m/s 
a 0 — 10 -6 m/s 2 
fa = 10- 6 N 

Another choice that is commonly adopted in instrumentation, e.g. spectrum 
analysers, is to take 1 V or the input voltage range (IVR) as the reference value. 
Vibration sensors, in fact, feed to the analyser an output voltage which depends — 
through the sensitivity of the sensor — on the vibration level to be measured. The 
analyser, in turn, allows the operator to preset an input voltage range 
(IVR = ±0.5 V, ±1 V, ±5 V, ±10 V (for instance) which is sometimes used 
as a reference when displaying power spectra on a dB scale. So, the dB value of 
a voltage signal fed to the analyser is given by 

L(dBV) = 201og 10 

and in this case one always gets negative dB levels, because always y m < IVR. 
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From the definition, it is clear that decibels are not added and subtracted 
linearly. If we have two decibel levels L 1 and L 2 (obviously referred to the 
same reference value) that must be added, the total dB level, say L r , is not 
L\ + Li but 

L t = 101og 10 (10 L| '' 10 + 10 Ll/1 °) (1.21) 





(b) 



Fig. 1.6 (a) Addition and (b) subtraction of dB levels (dB levels calculated as in eq 
(1.19)). 
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since 



10 L '/ ,0 = - and 10^ iO =- 
and the total level is by definition 




dB djftanrai L2-L1 (LZ>L1) 

(a) 



Fiji. 1 -7b: Subtradiari Of dB Level* [d6 level* calculated as in eq.(1-2DJl 




dB difference between total and background level* 

(b) 

Fig. 1.7 (a) Addition and (b) subtraction of dB levels (dB levels calculated as in eq 
( 1 . 20 )). 
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Similar considerations apply in the case of subtractions of dB levels, the 
most common case being the subtraction of a background level L b from a 
total measured level L x , in order to obtain the actual level L a . Of course, 
there is now a minus sign in the parentheses of (1.21) and in this case L\~Lj 
and Li — Ly . 

The graphs in Fig. 1.6(a), (b) (dB levels calculated according to eq (1.19)) 
and Fig. 1.7(a), (b) (dB levels calculated according to eq (1.20)) are useful 
when adding or subtracting dB levels and there is no pocket calculator at 
hand. For example, we want to add the two levels 1^=60 dB and L 2 =63 dB. 
Referring to Fig. 1.6(a), the difference A L = Li — Li is 3 dB and this value 
gives on the ordinate axis a dB increment of 1,76 = 1,8 to be added to the 
higher level (63 dB). The result of the addition is then 64.8 dB. 

For an example of subtraction, suppose that we measure a total level of 
70 dB on a background level of 63 dB. Referring to Fig. 1.6(b), for the 
abscissa entry AL — 7 dB °ne obtains for the ordinate value a dB decrement 
= 0.97 ^ 1.0 to be subtracted from the total level. The actual value is then 
69.0 dB. 

A final word on the calculation of the average dB level of a series of N 
measurements: this is given by 



L avg — 10 log 10 f — io L, A° 



and not by the familiar 



-I V' 

Lai'g ~ — ^ L, 



( 1 . 22 ) 



1.4 Springs, dampers and masses 

Almost any physical system possessing elasticity and mass can vibrate. The 
simplest realistic model considers three basic discrete elements: a spring, a 
damper and a mass. Mathematically, they relate applied forces to 
displacement, velocity and acceleration, respectively. Let us consider them 
more closely. 

The restoring force that acts when a system is slightly displaced from 
equilibrium is supplied by internal elastic forces that tend to bring the system 
back to the original position. In solids, these forces are the macroscopic 
manifestation of short-range microscopic forces that arise when atoms or 
molecules are displaced from the equilibrium position they occupy in some 
organized molecular structure. 

These phenomena are covered by the theory of elasticity. In addition, in 
complex structures, macroscopic restoring forces occur when different 
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structural elements meet. In general, it is not easy to predict the performance 
of a such a joint, even under laboratory conditions, because the stiffness is 
dependent on very local distortions. Here probably lies the main reason of 
the fallibility of predictions of stiffnesses in many cases. 

This may come as a surprise to many engineering graduates and it is due 
to the mistaken belief that ‘everything’ can be calculated with a high-speed 
computer, using finite-element techniques with very fine resolution. As it 
turns out, it is frequently cheaper to test than to predict. 

The easiest way to model such characteristics relating forces to 
displacements is by means of a simple stiff element: the linear massless spring, 
shown in Fig. 1.8. The assumption of zero mass assures that a force F acting 
on one end must be balanced by a force -F on the other end so that, in such 
conditions, the spring undergoes an elongation equal to the difference between 
the displacements x 2 and x t of the endpoints. For small elongations it is 
convenient and also correct to assume a linear relation, i.e. 

F ~ k(x2 — JCi ) (1-23) 

where k is a constant (the spring stiffness, with units N/m) that represents 
the force required to produce a unit displacement in the specified direction. 

Often, one end of the spring is considered to be fixed, the displacement of 
the other end is simply labelled x and — since F is an elastic restoring force — 
eq (1.23) is written F--kx, with the minus sign indicating that the force 
opposes the displacement. The reciprocal of stiffness, i.e. Ilk, is frequently 
used and is called flexibility or compliance. 

In real-world systems, energy is always dissipated by some means. In the 
cases concerning us the energies of interest are the kinetic energy of motion 
and the potential strain energy due to elasticity. Friction mechanisms of 
different kinds are often considered to represent this phenomenon of energy 
‘damping’ that ultimately results in production of heat and is probably the 
main uncertainty in most engineering problems. In fact, a word of caution 
is necessary: any claim that damping in structures can be predicted with 
accuracy should be treated with scepticism. 

On a first approach, we will see that damping can be neglected (undamped 
systems) without losing the physical insight to the problem at hand; 
nevertheless, it must be kept in mind that this is a representation of an unreal 
situation or an approximation that can be accepted only in certain 
circumstances. 



*1 *2 




Fig. 1.8 Ideal massless spring. 
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The simplest model of damping mechanism is the massless viscous damper 
that relates forces to velocities. An example can be a piston fitting loosely in 
a cylinder filled with oil so that the oil can flow around the piston as it 
moves inside the cylinder. The graphical symbol usually adopted is the dashpot 
shown in Fig. 1.9. Once again, a force F on one end must be balanced by a 
force -F on the other end and linearity is assumed, i.e. 



F — c(Xj — Jti ) (1-24) 

where the dot indicates the time derivative and c is the coefficient of viscous 
damping, with units N s/m. 

If one end is fixed and the velocity of the other end is labelled x, eq (1.24) 
is written F = —cx, with the minus sign indicating that the damping force 
resists an increase in velocity. 

The quantity that relates forces to accelerations is the mass and the 
fundamental relation is Newton’s second law, which can be written, with 
respect to an inertial frame of reference, 



F = mx (1-25) 

In the SI system of units the mass is measured in kilograms and eq (1.25) 
defines the unit of force (see also Section 1.1). The mass represents the inertia 
properties of physical bodies that, under the action of a given applied force 
F, are set in motion with an acceleration that is inversely proportional to 
their mass. 

Finally, it is interesting to consider a few examples of ‘equivalent springs’ 
k eq , meaning by this term the replacement of one or more combination of 
stiff elements with a single spring that takes into account and represents — 
for the problem at hand — the stiffness of such combination. 

Two springs can be connected in series or in parallel, as shown in Fig. 
1.10(a) and (b). (However, for a different approach see Section 13.7.3) In 
the first case we have 



F = kl(X2 — Xi) 
F = kiXi 



(1.26) 



For an equivalent spring of stiffness k eq connected between the fixed point 
and x 2 we would have F ~ k tq X 2 ■ It is then easy to show that it must be 



1 _ 1 1 _*! +k 2 



^4 



k-jkz 



-F 



-> 



(1.27) 



Fig. 1.9 Ideal massless dashpot. 
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and extending to the case of N springs connected in series we obtain 



k 



«<? — 




(1.28) 



For the connection in parallel of two springs (Fig. 1.10 (b)) the force F 
divides into F l and F 2 (with the obvious condition) F = Fj +Fi and the 
relations F\ = k-\x, F 2 =kix. The equivalent spring obeys F = k eq x, from 
which it follows that k eq = k\ + k±. The generalization to N springs in parallel 
is straightforward and we get 

N 

k eq =Y^ k > (1.29) 

; = i 



Considering the real spring of Fig. 1.11(a), under the action of an axial 
tension or compression its stiffness constant is 



Gd* 
k = 7 

64nR- 



(1.30) 



where G (in N/m 2 ) is the modulus of elasticity in shear of the material of 
which the spring is made, n is the number of turns and the other symbols are 
shown in Fig. 1.11a. Different values are obtained for torsion or bending 
actions and Young’s modulus, in these latter cases, also plays a part. 






i 



i 



I VVW\A WWW 

kj k 2 



(a) 



*2 



| WWW 

I 

p WWW 

/ 



-> F 



(b) 



X 



Fig. 1.10 Springs connected in: (a) series; (b) parallel. 
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Other examples are shown in: 

• Fig. 1 .1 1 (b): cantilever with fixed-free boundary conditions and force F 
applied at the free end, k is the local stiffness at the point of application 
of the force, 

• Fig. 1.11(c): fixed-fixed bar of length L with transverse localized force 
F at LI 2; k is the local stiffness at the point of application of the force; 

• Fig. 1.11(d): fixed-fixed bar of length L with uniform transverse load w 
(N/m), k is the local stiffness at the point of maximum displacement; 

• Fig. 1.11(e): bar simply supported at both ends with force F at LI 2, k is 
the local stiffness at the point of application of the force. 

The result of Fig. 1.11(b) comes from the fact that we know from basic 

theory of elasticity that the vertical displacement x of the cantilever free end 

under the action of F is given by 

_FL 3 
~ 3 El 



so that 

_F 3 El 

= x~U 

Similar considerations apply to the other examples, E being Young’s 
modulus in N/m 2 and I being the cross-sectional moment of inertia in m 4 . 

1.5 Summary and comments 

Chapter 1 introduces and reviews some basic notions in engineering vibrations 
with which the reader should already have some familiarity. At the very 
beginning, it is important to point out that the modelling process is always 
the first step to be taken in every approach to problems in science and 
engineering. The model must be able to reproduce, within an acceptable 
degree of accuracy, the essential features of the physical system under 
investigation. The definitions of the terms ‘acceptable degree of accuracy’ 
and ‘essential features’ depend on the particular problem at hand and, in 
general, should be based on decisions that specify: the results that are needed, 
the error one is willing to accept, the money cost involved, how and what to 
measure if measurements must be taken. When the budget is fixed and time 
deadlines are short, a compromise between cost and attainable level of 
accuracy must be made and agreed upon. 



Copyright © 2003 Taylor & Francis Group LLC 



More specifically, a few classifications are given which may help in setting 
the guidelines of the type of modelling schemes that can be adopted for a 
wide class of problems. 

First and above all is the distinction between linear and nonlinear analysis. 
Linearity or nonlinearity are not intrinsic properties of the system under 
study, but different behaviours of mechanical and structural systems under 
different conditions. Small amplitudes of motion, in general, are the range 
where linearity holds and the cornerstone of linearity is the principle of 
superposition. 

Second, discrete and continuous systems can be distinguished, or, 
equivalently, finite or infinite-degree-of-freedom systems. Continuous 
distributed parameters are often substituted by discrete localized parameters 
to deal with ordinary differential equations rather than with partial 
differential equations and perform the calculations via computer. 

Third, signals encountered in the field of linear vibrations can be classified 
as deterministic or random: an analytical form can be written for the former, 
while statistical methods must be adopted for the latter. The type of signals 
encountered in a particular problem often dictates the method of analysis. 

The other sections of the chapter introduce some basic definitions and 
methods that will be used throughout the text. A few examples are simple 
sinusoidal motion, its complex (phasor) representation and the decibel scale. 
The phenomenon of beats is then considered, for its own intrinsic interest 
and as an application of phasors, and, finally, an examination of the 
parameters that make systems vibrate or prevent them from vibrating — 
namely mass, stiffness and damping — is made, together with their simplest 
schematic representations. 
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2 



Mathematical preliminaries 



2.1 Introduction 

The purpose of this chapter is twofold: ( 1 ) to introduce and discuss a number 
of concepts and fundamental results of a mathematical nature which will be 
used whenever necessary in the course of our analysis of linear vibrations, 
and (2) to provide the reader with some notions which are important for a 
more advanced and more mathematically oriented approach to the subject 
matter of this text. In this light, some sections of this chapter can be skipped 
in a first reading and considered only after having read the chapters that 
follow, in particular Chapters 6-9. 

It is important to point out that not all the needed results will be considered 
in this chapter. More specifically, matrix analysis is considered separately in 
Appendix A, while the whole of Chapter 1 1 is dedicated to the discussion of 
some basic concepts of probability and statistics which, in turn, serve the purpose 
of introducing the subjects of stochastic processes and random vibrations (Chapter 
12). Also, when short mathematical remarks do not significantly interfere with 
the main line of reasoning of the subject being considered, brief digressions are 
made whenever needed in the course of the text. 

So, without claim of completeness and/or extreme mathematical rigour, 
this chapter is intended mainly for reference purposes. We simply hope that 
it can be profitably used by readers of this and/or other books on the specific 
field of engineering vibrations and related technical subjects. 

2.2 Fourier series and Fourier transforms 

In general terms, Fourier analysis is a mathematical technique that deals 
with two problems: 

1. the addition of several sinusoidal oscillations to form a resultant 
(harmonic synthesis); 

2. the reverse problem, i.e. given a resultant, the problem of finding the 
sinusoidal oscillations from which it was formed (harmonic analysis). 
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As a trivial opening example, it is evident that adding two harmonic 
oscillations of the same frequency 

JCi = Ai sin wot + Bi costv'o? 

( 2 . 1 ) 

,X2 = A 2 sinwQ? + Bi cosuot 

leads to a harmonic oscillation xq + xi of the same frequency with a sine 
amplitude Ai +Az and a cosine amplitude B] +Bn. If the two oscillations 
have different frequencies, the resultant will not be harmonic. However, if 
the frequencies being summed are all multiples of some fundamental frequency 
coi, then the resultant oscillation will repeat itself after a time T = 2tt/cji and 
we say that it is periodic with a period of T seconds. For example, the function 

x(t ) — cos (3i) + 1.4sin(5f) + 2.1 cos (6t) (2-2) 

shown in Fig. 2.1, repeats itself with period T — 2 tt/1 — 2tt seconds, so that 
x(t) = x{t + T) (in Fig. 2.1, note that the time axis extends from j — 0 to 
t = 6ir^ 18.85). 

If, on the other hand, the frequencies being summed have no common 
factor, the resultant waveform is not periodic and never repeats itself. As an 
example, the function 

x(f) = cos(t) + l^sinfv^i) + 1.3cos(\/l7f) (2-3) 

is shown in Fig. 2.2 from f=0 to t=50 seconds. 




Fig. 2.1 Periodic function. 
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Fig. 2.2 Nonperiodic function. 

In short, the process of harmonic synthesis can be expressed 
mathematically as 

N 

x(r) = ^(A„ costu„r + B„ sin(u n t) (2.4) 

nr- 1 

where the quantities A„, and B„ are, respectively, the cosine and sine 
amplitudes at the frequency co n . When a fundamental frequency tOj exists, 
then — fUO) and we call a>„ the nth harmonic frequency. 

2.2.1 Periodic functions: Fourier series 

As stated before, harmonic analysis is, in essence, the reverse process of 
harmonic synthesis. J.B. Fourier’s real achievement, however, consists of the 
fact that he showed how to solve the problem by dealing with an infinite 
number of frequencies, that is, for example, with an expression of the type 

_OQ_ 

x(t) = ^[A„ cos(nwit) + fl„ sin(«wit)] ( 2 . 5 ) 

M ~ 1 



where it is understood that the meaning of eq (2.5) is that the value of x(t) 
at any instant can be obtained by finding the partial sum S N of the first N 
harmonics and defining x(t) as the limiting value of S N when N tends to 
infinity. As an example, it is left to the reader to show that eq (2.5) with the 
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coefficients A„= 0, and B„=l/n is the well-known periodic ‘sawtooth’ oscillation 
(see Fig. 2.3 representing the partial sum S 6 with Wj = 1 rad/s). 

Incidentally, two short comments can be made at this point. 

First, the example of Fig. 2.3 gives us the opportunity to note that a 
function which executes very rapid changes must require that high-frequency 
components have appreciable amplitudes: more generally it can be shown 
that a function with discontinuous jumps (the ‘sawtooth’, for example) will 
have A and B coefficients whose general trend is proportional to n- 1 . By 
contrast, any continuous function that has jumps in its first derivative (for 
example, the triangular wave of Fig. 1.5) will have coefficients that behave 
asymptotically as n~ 2 . 

Second, although our notation reflects this particular situation, the term 
‘oscillation’ does not necessarily imply that we have to deal with time-varying 
physical quantities: for example, time t could be replaced by a space variable, 
say z, so that the frequency co would then be replaced by a ‘spatial frequency’ 
(the so-called wavenumber, with units of rad/m and usually denoted by k or k), 
meaning that x(z) has a value dependent on position. Moreover, in this case, the 
quantity 2n/k (with the units of a length) is the wavelength A of the oscillation. 

So, returning to the main discussion, we can say that a periodic oscillation 
x(t) whose fluctuations are not too pathological (in a sense that will be 
clarified in the following discussion) can be written as 





(2.6a) 
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Fig. 2.3 Harmonic synthesis of sawtooth oscillation (first six terms). 
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and the Fourier coefficients A n , and B„ are given by 

2 

A„ = — x(t)cos(JiW[t)df 

T J!l 

2 C (2.6b) 

B„= — ;e(f)sin(wwit)d£ 

T Jti 



where 2-jt/uq — T is the period of the oscillation and the limits of integration 
£ b t 2 can be chosen at will provided that ti~ t\=T (the most frequent choices 
being obviously t\~Q,ti = T or t\ — — T /2, ti — T (2). Note that the ‘static’ 
term (1/2 )A 0 has been included to allow the possibility that x(t) oscillates 
about some value different from zero; furthermore, we write this term as (1/ 
2 )A 0 only for a matter of convenience. By so doing, in fact, the first of eqs 
(2.6b) applies correctly even for this term and reads 

Ao 1 f T 

— = - x[t)dt (2.6c) 

2 T Jo 

where we recognize the expression on the r.h.s. as the average value of our 
periodic oscillation. 

As noted in Chapter 1, it often happens that complex notation can provide 
a very useful tool for dealing with harmonically varying quantities. In this 
light, it is not difficult to see that, by virtue of Euler’s equations (1.7), the 
term A cos ujt + B sin wt can also be written as 



C ' e‘~' ■ C e ! ~‘ (2-7a) 

where the complex amplitudes C + and C” are given by 
. 1 

C + =-(A-iB) 

2 

1 (2.7b) 
C” = - (A + j'B) 

2 



and no convention of taking the real or imaginary part of eq (2.7a) is 
implied because C + and C“ combine to give a real resultant. Then, eq (2.6a) 
can also be written as 



*(£) = * A 0 + f] (C:e^ + C- n e—“) 

2 H-l 



DC 




n = - dc 



( 2 . 8 ) 
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where in the last expression the complex coefficients C„ are such that C„ = C%, 



1 

C„ =- (A u - iB n ) 
2 

C n =- A 0 

2 

1 



C~„ — — (A„ + iB„) 
2 



(2.9) 



and eqs (2.6b) translate into the single formula 




-T 

o 



xUie-^'dt 



( 2 . 10 ) 



Now, if we become a bit more involved with mathematical details we can 
turn our attention to some issues of interest. The first issue is how we obtain 
eqs (2.6b) (or eq. (2.10)). Assuming that our periodic function x(t) can be 
written in the form 



x[t)= 

n — - -X) 

let us multiply both sides of this equation by and integrate over one 

period, to get 

tT »c ,t 

x{t)e~‘ m * l dt= V (2.11) 

JO n — - oc JO 

Noting that 

\ T e ! ’ t ^ t e-' m ^ t dt=Td mn 

Jo 

where 8„,„ is the Kronecker delta defined by 







m^n 
m = n 



then 



■T 

x{t)e~ i ” l ^ 1 dt=TC n 6 mn 



Jo 






( 2 . 12 ) 



which is precisely eq (2.10). 
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The second issue we consider has to do with the conditions under which 
the series (2.6a) (or (2.8)) converges to the function x(t). As a matter of fact, 
three basic assumptions have been tacitly made in the foregoing discussion: 
(1) the expansion (2.6a) is possible, (2) the function under the 

integral in eq (2. 1 1 ) is integrable over one period and (3 ) the order of integration 
and summation can be reversed (r.h.s. of eq (2.11)) when calculating the Fourier 
coefficients (eqs (2.6b) or (2.10)). Various sets of conditions have been 
discovered which ensure that these assumptions are justified, and the Dirichlet 
theorem that follows expresses one of these possibilities: 

Dirichlet theorem. If the periodic function x(t) is single valued, has a finite 
number of maxima and minima and a finite number of discontinuities, and 
if Jq | x(r) | dt is finite, then the Fourier series (2.6a) converges to x(t) at all 
the points where x(t) is continuous. At jumps (discontinuities) the Fourier 
series converges to the midpoint of the jump. Moreover, if x(t) is complex (a 
case of little interest for our purposes), the conditions apply to its real and 
imaginary parts separately. 



Two things are worthy of note at this point. First, the usefulness of this 
theorem lies mainly in the fact that we do not need to test the convergence 
of the Fourier series. We just need to check the function to be expanded, and 
if the Dirichlet conditions are satisfied, then the series (when we get it) will 
converge to the function x(t) as stated. Second, the Dirichlet conditions are 
sufficient for a periodic function to have a Fourier series representation, but 
not necessary. In other words, a given function may fail to satisfy the Dirichlet 
conditions but it may still be expandable in a Fourier series. 

The third and last issue we consider is the relation between the mean 
squared value of x(t), i.e. 




r t 

0 



x z {t)dt 



(2.13) 



and the coefficients of its Fourier series. The result we will obtain is called 
Parseval’s theorem and is very important in many practical applications where 
energy and/or power are involved. If, for example, we use the expansion in 
exponential terms (eq (2.8)) it is not difficult to see that 




dt 



(2.14a) 
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which shows that each Fourier component of x(t), independently of the other 
Fourier components, makes its own separate contribution to the mean squared 
value. In other words, no cross (or ‘interference’) terms of the form 
C„C* ; (» 5 ^ m) appear in the expression of the mean squared value of x(t). If, 
on the other hand, we use the series expansion in sine and cosine terms, it is 
immediate to determine that Parseval’s theorem reads 

(x 2 (t))=^- + ^ ( A l+ B2 J (2.14b) 

4 2 „ = i 



2 . 2.2 Nonperiodic functions: Fourier transforms 

Nonperiodic functions cannot be represented in terms of a fundamental 
component plus a sequence of harmonics, and the series of harmonic terms 
must be generalized to an integral over all values of frequency. This can be 
done by means of the Fourier transform, which, together with the Laplace 
transform, is probably the most widely adopted integral transform in many 
branches of physics and engineering. 

Now, given a sufficiently well-behaved nonperiodic function f(t) (i.e. the 
counterpart of the function x(t) of Section 2.2.1) the generalization of the 
Fourier series to a continuous range of frequencies can be written as 

f{t)= f F(cj)e" ! <*w (2.15) 



provided that the integral exists. Then, the counterpart of eq (2.10) becomes 

1 r 00 

F(w)=— f(t)e~^dl (2.16) 

2tt J -oo 

where the function F(co) is called the Fourier transform of f(t) and, conversely, 
f(t) is called the inverse Fourier transform of F(co). Also, the two functions 
f(t) and F(co) are called a Fourier transform pair and it can be said that eq 
(2.15) addresses the problem of Fourier synthesis, while eq (2.16) addresses 
the problem of Fourier analysis. 

A set of conditions for the validity of eqs.(2.15) and (2.16) are given by 
the Fourier Integral theorem, which can be stated as follows: 

Fourier integral theorem. If a function f(t) satisfies the Dirichlet conditions 
on every finite interval and if j fit) \ dt < do, then eqs (2.15) and (2.16) 
are correct. In other words, if we calculate F(a>) as shown in eq (2.16) and 
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substitute the result in eq (2.15) this integral gives the value of f(t) where f(t) 
is continuous, while at jumps, say at t-t 0 , it gives the value of the midpoint 
of the jump, i.e. the value [(/fr^ ) +f(t g )]/2. 

Without giving a rigorous mathematical proof (which can be found in 
many excellent mathematics texts), let us try to justify the formulas above. 
Starting from the Fourier series of the preceding section, it may seem 
reasonable to represent a nonperiodic function f(t ) by letting the interval of 
periodicity T— »oo. So, let us rewrite eq (2.8) as 

30 PC 

m= £ c ^ ,! = £ ( 2 .i 7 ) 

- 00 — 30 

where w K = «W] and we know that the fundamental frequency is related 
to the period T by W] = 2ir/T. In this light, if we define Aw = w„ + i — w„, it 
is seen immediately that Aw = 2? t/T so that eq (2.10) becomes 

Aw fT/2 

C n = f(t)e ,Wa 'dt (2.18) 

J-t/2 

and can be substituted in eq (2.17) to give 



fit) 





f(a)e 



e^'Aw 



1 

2tt 



£ 



/ j-T/2 \ 1 00 

! °^da> )Aw = — g{w„) Aw (2.19a) 

y J - T /2 J 2tt “ 



where we define (note that a is a dummy variable of integration) 



fT/2 

gK) = f{a)e Mn{t ~ a} da 
i-T/2 



(2.19b) 



Now, if we let T —■ » oo, then Aw — » 0 and, in the limit 



1 °° 1 
— £ gba n )Aw— » — 
2 * 2it J 



g{uj)du> 



gi w«) 



fia^-^da 
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so that eq (2.19a) becomes 






J-= 



r x l 

g(u)duj = — 

-oo 2 tt 



f(a)e^ {t o] dadu} 




'da e^du = 



F(oj}e’^dLj 



( 2 . 20 ) 



where we have defined 



m 



1 r 00 

= — f{a)e~ 

2.7T 1— oc 



'da 



( 2 . 21 ) 



Equations (2.20) and (2.21) are the same as (2.15) and (2.16). Equation 
(2.20) also justifies the different notations frequently encountered in various 
texts where, for example, the multiplying factor 1/(2jc) is attached to the 
inverse transform. The difference is irrelevant from a mathematical point of 
view, but care must be exercised in practical situations when using tables of 
Fourier transform pairs. Some commonly encountered forms (other than 
our definition given in eqs (2.15) and (2.16)) are as follows: 



1 

m =— 

2i r 



fM = 



F(u)e ! “'*du 



f{t)e~‘“ ! dt 



(2.22a) 



m = 



J 



F^e'^dui 



F(w)=— = [ f{t)e 
y/2n i-x 



*dt 



(2.22b) 



f(t)= p F(v)e aimt dv 
J--x 

F(i>)= [ f{t)e ,2lTi,l dt 

J—act 



(2.22c) 



where in eqs (2.22c) the ordinary frequency v (in Hz) is used rather than the 
angular frequency w = 2w (in rad/s). 
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Example 2.1 Given a rectangular pulse in the form 



/(*) 



1/2 -1 < r < 1 

0 otherwise 



(2.23) 



(a so-called ‘boxcar’ function) we want to investigate its frequency content. 
From eq (2.16) we have 



1 

EM = — 

2tt 



fi 

J-i 



1 

= ■ (e 

4 mui 





4mu) 



sinw 
2i rw 




(2.24a) 



and the graph of F(co) is sketched in Fig. 2.4. 

Furthermore, if we want the Fourier transform in terms of ordinary 
frequency, from the condition F(uj)duj = F{v)dv and the fact that du — 2 xdv 
we get 



F{v) = 2 ttE(w) 



sin(2-?n/) 

2x1^ 



(2.24b) 



Note that in this particular example F(co) is a real function (and this is 
because the real function f(t) is even). In general, this is not so, as the reader 
can verify, for instance, by transforming the simple exponential function 




Fig. 2.4 Fourier transform of even ‘boxcar’ function. 
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f{t) = exp(— t), which we assume to be zero for t< 0. The result is 

1 - (U> 

fM = 

2tt(1 + lj 2 j 

Alternatively, the reader can consider the Fourier transform of the pulse 
function f(t)= 1/2 for 0<t<2 and zero otherwise (which is just a time-shifted 
version of (2.23)) and note that the result can be written as 

e -iw 

F(iu) = sincj 

2jtW 

with exactly the same magnitude as the expression (2.24a) but different phase 
(because the phase depends on the arbitrary definition of time zero). 

Before considering some properties of Fourier transforms, a word on notation 
is necessary: in what follows — and whenever convenient in the course of 
future chapters — the symbol F(w) = will indicate that F(co) is the 

Fourier transform of f(t). Conversely, the symbol f{t) — S' -1 (FM) will indicate 
that f(t) is the inverse Fourier transform of F(a>). With this notation, the first 
easily verified property is linearity, that is 

9{af(t) + bg(t)}=a9{f(t)} +&${*(*)} (2.25) 

where a and b are two constants and f(t) and g(t) are two Fourier- 
transformable functions. 

A second important property has to do with the Fourier transform of the 
function df(t)/dt. Integrating by parts and taking into account the fact that — for 
the transform to exist — necessarily we must have f(t) —> * 0 as t * ±co, we get 






{ 




1 f 00 


df 

— - e~‘^dt 


2?r J-oo 


dt 


i 


isjj 


— [fe 


“* 1 “ + — 
J — CC 1 „ 


27T 


2tt 







(2.26a) 



which is just a special case of the relation 






{ 



d"f(t) 

dt n 






(2.26b) 



Clearly, eq (2.26b) is true only if d”f jdt n is a Fourier-transformable function 
in its own right. By similar arguments, if f(t) is absolutely integrable then the 
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function /(t) = J) fi'yjd'j is continuous. If it is also absolutely integrable, 
then 



W)} = — 

2tv 



Ie-^dt 



1 — T fe~ i *dt = -9{f(t)} 

Inivj J oe 



2k iuj 



iu 



(2.27) 



Consider now two Fourier-transformable functions f(t) and g(t) with Fourier 
transforms F(co) and G(co), respectively. The function w(t) defined by 



Mt) = f(t) *g(t) = 



f(t - T)g{T)dft 



(2.28) 



is called the convolution of f(t) and g(t) and it has considerable importance 
when calculating the time response of linear systems. Provided that 
J™oe I | dt < oo we can take the transform of w(t) and use Fubini’s 

theorem (see any book on calculus) to get 



1 

2ft 



f(t — r}g(T)e drdt 



1 

27T 



■ oc 




r sc 1 






1 

1 

t 

£ 


OO 


* 


f -oc J 



dr 



Then, by making the change of variable <p = t — r, dip = dt in the t integral, 
it follows that 



W(w) 



■ j : 






1 

2i r 



f{p)e ‘^dtp 



dr 



= F(w) g(T)e~^dT = 2nF(u)G(u) 



(2.29a) 



which expresses the convolution theorem: in words, the transform of the 
convolution f(t) * g(t) is 2n times the product of the transforms of f(t) and 
g(t). Again, the position of the 2n factor depends on the definition of Fourier 
transform that we use (eqs (2.22a-c)). By the same token, it is left to the 
reader to show that 

^- 1 {TM*G( w ))=/ r (%(t) (2.29b) 
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The Parseval formula for Fourier transforms reads 



| f{t) 1 2 dt = 2-tt 



| F(w)| 2 duj 



and can be derived as follows. Consider the integral 



_ 1 
2tt J 



f{t)F*{u)e~ ,wt dtdu 



(2.30) 



Assuming that the order of integration is immaterial (again Fubini’s theorem), 
we can write it as 



/ = 



■oo ^ poo 

F*M — f{t)e~ iw> dt 
J-oc 2.7T J— oc 






= J™ F*{Lj}Fdui= | |F(w)| 2 dw 



(2.31) 



or, noting that f(t) = & 1 {F(uj)} and hence f *(f) = F*(oj)e 



as 



1 = 



F *(u})e~ Llt du 



— r ^(^*(0*=^- r i^wr 

2tt J-w 2tt J-oc 



dt 



(2.32) 



so that eqs (2.31) and (2.32) together lead exactly to the Parseval theorem 
expressed by eq (2.30). Like its counterpart for periodic functions (eqs (2.14a) 
or (2.14b)), eq (2.30) is very important in practical applications when dealing 
with the squares of physical time-varying signals, i.e. quantities that have to 
do with the energy and/or the power associated with the signal itself. 



2.2.3 The bandwidth theorem ( uncertainty principle) 

The idea of the uncertainty principle comes from the development of quantum 
mechanics during the 1920s. The principle, in its most famous form given 
by the physicist W. Heisenberg, reads AxAp > h and, broadly speaking, means 
that the product of the uncertainties in position (Ax) and momentum (A p) 
for a particle is always greater than Planck’s constant 6—1.05 x 10 _34 Js. 
More generally, the principle states that any phenomenon described by a 
pair of conjugate (or complementary) variables must obey some form of 
uncertainty relation. 
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For our purposes the conjugate variables are time t and frequency v and 
the uncertainty principle in the field of signal analysis relates the duration of 
a transient function (a pulse) to the range of frequencies that are present in 
its transform. So, if we call At the duration of the transient f(t ) and Av the 
range of frequencies spanned by F(v), then 

AtAi/^rl (2.33) 



which is the so-called bandwidth theorem. The approximation sign in eq 
(2.33) means that the product lies usually in the range 0.5-3 for most simple 
transients (try, for instance, the boxcar function of Example 2.1, where Av 
can be approximately taken as the width at the basis of the peak centred at 
v=0), but its precise value is not really important for the essence of the 
argument. Stated simply, eq (2.33) shows that the two members of a Fourier 
transform pair — each one in its appropriate domain — cannot both be of short 
duration. The implications of this fact pervade the whole subject of signal 
analysis and have important consequences in both theory and practice. 

In the course of this text, we will encounter many situations that agree 
with eq (2.33): for example, a lightly damped structure in free vibration will 
oscillate for a long time (large At) at its natural frequency m„, so that the 
spectrum of the recorded signal will show a very narrow peak at ui = 
(small Am); by contrast, if we want to excite many modes of a structure in 
a broad band of frequencies (large Am), we can do so by a sudden blow of 
very short duration (small At) or by means of a random, highly erratic time 
signal, and so on. So, the essence of this short section is that the uncertainty 
principle represents an inescapable law of nature (or, at least, of the way in 
which we perceive the phenomena of nature) with which we must come to 
terms either by looking for trade-offs or by using different methods to extract 
all the desired information from a given signal. 

Finally, if we observe that for slowly decaying signals it may not be so 
easy to define the quantities At and Av, we can give the following definitions 
in order to make these concepts more precise: 



(At) 2 



‘00 

t z \f(t)\ 2 dt 

J— 00 



\m\ 2 dt 

— 00 



( A ^) 2 



v 2 | F{v) \ 2 dv 

J -00 



‘00 

I F(v) 1 2 dv 

J -OG 



(2.34) 
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These integrals may not be easy to calculate, but they give a prescription 
applicable to a wide variety of signals: with these definitions it can be shown 
that the bandwidth theorem becomes AtAi/> 1/(4^}. 



2.3 Laplace transforms 



From a general mathematical viewpoint, the Fourier transform is just a special 
case of a class of linear integral transforms which can be formally written as 






■b 



K(t , u)f{t)dt 



(2.35) 



where the function K(t, u) is called the kernel of the transform and T{f(t)} is 
the integral transform of f(t) with respect to the kernel K(t,u). The various 
transforms are characterized by their kernels and their limits a and b so that, 
for example, for the Fourier transform K(t, w) = (l/(2ir)}e ~ mt and U — — oc, 
b = o o. Other types of integral transforms are, to name just a few, the Laplace 
transform, the Flankel transform and the Mellin transform. 

Mathematically, their utility lies in the fact that, with certain kernels, 
ordinary differential equations yield algebraic expressions in the transformed 
variable (partial differential equations also lead to more tractable problems) 
and the solution in the transformed space is found more easily. Clearly, the 
problem of finding the solution in the original variable remains and the 
calculation of the inverse transformation is often the most difficult part of 
the whole procedure. 

In its own right — together with the Fourier transform — the Laplace 
transform is the most popular integral transform. It is defined as 



f{s) = <£{/(()}= [ f{t)e~«dt (2.36) 

Jo 

where s is a complex quantity. Equation (2.36) is particularly useful when 
we are interested in functions whose Fourier transform does not exist, i.e. 
for example the function f(t) = e b> , where b> 0. In many cases the trouble at 
£ — > 00 can be fixed by multiplication with a factor e~ ct if c is real and larger 
than some minimum value a (which, in turn, depends on the function to be 
transformed: for f(t)-e bt , for example, we must have ob), while the ‘bad’ 
behaviour at £ — » — $o can be taken care of by noting that our interest often 
lies in functions f(t) that are zero for t< 0. In these circumstances, the function 
g(£) = f{t)e~ ct does have a Fourier transform and we get 

fX f 50 - 

f[t)e~ u e~ Llt dt = f(t)e~ {c + lu;]t dt = F(c + iu) 

,o Jo 
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Introducing the complex variable s — c + iw, we obtain exactly the integral 
of eq (2.36) which, in turn, exists only in the right half of the s-plane where 
c = Re(s) > a (a is the minimum value for c mentioned above). 

Therefore, broadly speaking, we can say that the Laplace kernel is a 
damped version of the Fourier kernel or, more properly, that the Fourier 
kernel is an undamped version of the Laplace kernel. The inverse transform, 
in turn, can be written as 

1 f 00 

f{t)e~ a = — F(c + (2.37) 

27T J —oo 



Now, since s = c + iuf and hence du> = d$/i, eq (2.37) becomes 



f(t)e- a = — 
2m 



'C + iac | 

F(s)e {s - C}t ds = ~ e~ c ‘ 

c-iso 2m 



F(s)e s ‘ds 



from which it follows that 

J[ f c + ioo 

X~ 1 {F[s))sf(t)=— F(s)e s, ds (2.38) 

2m T- - ix 



where the integral on the r.h.s. is known as Laplace inversion (or Bromwich) 
integral and is understood in the principal-value sense. It converges to f(t) 
where f(t) is continuous, while at jumps it converges to the midpoint of the 
jump; in particular, for t - 0 the integral converges to (1/2)/" (0 + ). 

Note that we have made the factor l/(2n) appear in the inverse transform 
and not in the forward transform; this is merely a matter of convenience 
due to the fact that the Laplace transform is almost universally defined as 
in eq (2.36). 

Although in practical situations one generally refers to tables of Laplace 
transform pairs (e.g. Erdelyi [1 J), the integral of eq (2.38) can be evaluated 
as a complex (contour) integral by resorting to the theorem of residues and 
the following understanding: 

1. For £>0 the integration is calculated along the straight vertical line 
Refs)=c, where c is large enough so that all the poles of the function 
F(s)e st lie to the left of this line. The contour is closed, forming a large 
semicircle around the complex plane to the left of the line. 

2. For £<0, since f(t) must be zero, the contour is closed, forming a semicircle 
that extends to the right of the line (no poles in the contour). More 
details on complex integration can be found, for example, in Mathews 
and Walker [2] or Sidorov et al. [3]. 
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Now, if we turn to some properties of Laplace transforms, it is not difficult 
to show that, first of all, the Laplace transform is a linear transformation 
and that the first derivative of f(t) transforms into 






df{t) 

dt 



= sF{s)-m 



Also, a double integration by parts leads to 






{ 




s 2 F(s)-sf{ 0)~ 



df{t) 

dt 



i 



o 



or, if we write 

f w it) = df{t)jdt ., f <">(*) = d*m/dt* 

then, more generally 





n— 1 

V sT^-'-^O) 


l. / 


_ ; = 0 



(2.39) 



(2.40) 



(2.41) 



where — in eqs (2.39), (2.40) and (2.41) — it must be noted that 0 really means 
0 + , the limit as zero is approached from the positive side. 

If f(t) is piecewise continuous, then the function I(t) = jjj f(a}da is 
continuous, I( 0)=0 and we get 



£e{J(r)] =£< 



f(cx)da 






(2.42) 



Furthermore, given two functions f(t) and g(t), the convolution theorem 
for Laplace transforms reads 



* gW} = 2{/W}2fe(f)} = F(s)G(s) (2.43) 

where in the rightmost expression we have defined F(s) and G(s), the Laplace 
transforms of f(t) and g(t), respectively. 



Example 2.2. With reference to the problem of finding the inverse Laplace 
transform of a function it is of interest to note that one often has to deal with 
functions which can be written as the ratio of two polynomials in s, i.e. 



iW 



P(s) 

QU) 



(2.44a) 
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where Q(s) is of higher degree than P(s). If we suppose that Q(s) is a polynomial 
of order n with the n distinct roots qi, qz-, ■■■, q n , then the polynomial in the 
denominator can be written as Q(s) = (s — — qz) ... (s — q „ ) and the 

function F(s) as 

" Aj 

F i s ) = 2s (2.44b) 

; = l s 

where the coefficients A ; are obtained from 

Aj = lim [(s - ?/)F(sl] (2.44c) 

S— ' 



Since 







= exp(^f) 



(2.45) 



the inverse transform of F(s) is 

/‘(0 = S _1 JF(s)}= A, exp(qjt) 



(2.46) 



As a simple example, consider the function 






5+1 
s(s + 2) 



where the roots of the denominator are qi = 0, qi = — 2. From eq (2.44c) we 
get A\ = Ai = 1/2; therefore 

1 1 

F(s)=- + 

25 2(s + 2 ) 



and 



, 1 1 
fit) = - + -e 
2 2 



- 2 1 



As for the case of repeated roots and/or more difficult cases in general, 
the reader can find more details on partial-fraction theory in any basic text 
on calculus or algebra. 
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In the light of the foregoing discussion, we may observe that the main 
advantages of the Laplace transform over the Fourier transform are that: 

1. The Laplace integral converges for a large class of functions for which 
the Fourier integral is divergent. 

2. By virtue of eq (2.41), the Laplace transform provides automatically 
for the initial conditions at t- 0. This latter characteristic will be 
considered (and exploited) in the following examples and in future 
chapters. 



Example 2.3. For this example we consider the ordinary (homogeneous) 
differential equation 



d z f{t) 
dt 2 



+ a 2 f(t) = 0 



(2.47) 



where a is a constant and we are given the initial conditions f(0)=fo and 
df /dt | ( = 0 — f 0 . Taking the Laplace transform on both sides gives 

[s 2 F(s)-s/ r 0 -f o ']+j 2 F(s)=0 

so that the solution in the s-domain is easily obtained as 

m — S h-+Ji- 

s 2 +a 2 s 2 +a 2 



Then, from a table of Laplace transforms, we get % 1 [s/(s 2 + a 2 }} = cos at 
and 2? _1 {(s 2 +£ 2 ) -1 } = a~ l smat and hence 



fo ■ 



f{t ) = /ocos at H — sin^t 
a 



(2.48) 



which is precisely the result that we will obtain in eqs (4.8) and (4.9) by 
standard methods. 

On the other hand, if we consider the nonhomogeneous differential 
equation 



d 2 f(t) 
dt 2 



+ a 2 f(t) =g(t) 



(2.49) 



where, for example, the forcing function is g(f)=gcosutf and the initial 
conditions are still given by f( 0) = fo, df/dt |, _ 0 = f 0 , Laplace transformation 
of both sides leads, after a little manipulation, to 
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u 



f{s) = 



gs 



s f() 



■ + ■ 



(s- +a 2 ){s 2 + w 1 ) ($ 2 + tf 2 



because 5£{cosu;f} = s/{s 2 + uj 2 ), In order to return to the time domain, we 
already know the inverse transform of the last two terms, while for the first 
term we can use the convolution theorem and write 



5T 



s 2 + u> 2 } \ s 2 + a 2 



= a 



-l , 



1 



s 2 +u 2 



*2 



-i 



= [cos wf] * [a 'sin at] = - 



a Jo 



cos(iur)sin a{t 



so that the time domain solution of eq (2.49) is 

cos wr sin a(t — r)dr 



a . i " 



/(f) — fo cos at H — sin at + - 
a a 




- T)dr (2.50) 



(2.51) 



(see also eq (5.19) with the appropriate modifications). 

Partial differential equations can also be solved with the aid of Laplace 
transforms, the effect of Laplace transformation being the reduction of 
independent variables by one. Also, Laplace and Fourier transforms can be 
used together, as in the following example. 



Example 2.4. (Initial- value problem of an infinitely long flexible string). This 
problem will be considered in detail in Section 8.2; see this section for the 
meaning of the symbols. 

The equation of motion for the small oscillations of a vibrating string is 



d z y 1 d 2 y 
dx 2 c 2 dt 2 



(2.52) 



Let the initial conditions of the problem be specified by the two functions 
y(x, 0) = u{x) 



8y(x, t ) 



dt 



— w(x) 



2=0 



(2.53) 
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If we call Y(x, s) the Laplace transform of y(x, t) with respect to the time 
variable and transform eq (2.52), we get 

d 2 Y 1 , 

= — [s 2 Y — swfx:) — iv(x)] (2.54) 

dt 2 c 2 

Now, let us take the Fourier transform of eq (2.54) with respect to the 
space variable x and define 

tf(M = W*,s)} 

U(k) = '&{u(x}} 

W{k) = &{w{x)} 



where the variable k (the wavenumber) is the conjugate variable of x. We have 
s 2 s 1 

- * 2 *(M =- mk,s) - ~ u{k ) - - w(k> 

c 2 c 2 c 2 

from which it follows that 






sU(k) + W(k) 
s 2 +k 2 c 2 



(2.55) 



Taking the inverse Laplace transform of (2.55) gives 

, W 

X(k,t)=% 1 {(&(£, s)} = t/(£)cos(&:f) 4 sin(^ct) (2.56) 

kc 

and the final solution can be obtained by inverse Fourier transformation of 
eq (2.56). The inverse Fourier transform of the first term on the r.h.s. is 



&-'{U{k}cos(kct)} 



1 

2 



'OO 

U{k)[e ,kct + e- ihct ]e lkx dk 

J -QC 



1 

2 




U{k)e iHx + ct) dk 




U{k)e‘ k{x - Ct) dk 



1 

2 



[w(* + ct) + u(x — c£)] 



(2.57a) 
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while for the second term we may note that 



;-1 



W(k) 



sin(&ct) 1 



kc 






cos(&ca)da j 



where a is a dummy variable of integration. Hence 



S'M W(k) 



sin(&ci) \ 



kc 



-hi 



1 + 






(2.57b) 



and putting eqs (2.57a) and (2.57b) back together yields 

1 1 f* + <* 

y(x,t) = '9 = - M* + ct) + u{x - ct)\ -1 

2 2 c Ji - a 



(2.58) 



which is the same result we will obtain in eq (8.7a) by a different method. 

Incidentally, note that for a string of finite length L, we can proceed from 
eq (2.54) by expanding the functions Y(x, s), u(x) and w(x) in terms of 
Fourier series rather than taking their Fourier transforms. 



2.4 The Dirac delta function and related topics 

In physics and engineering it is often convenient to envisage some finite 
effect achieved in an arbitrarily small interval of time or space. For example, 
the cases with which we will have to deal in future chapters are impulsive 
forces which act for a very short time or localized actions applied at a given 
point on a structure. 

In this regard, the general attitude consists of using mathematics as a 
means for handling relations reasonably and efficiently, rather than as a 
free-standing discipline. Under the guide of mathematical reasonableness and, 
hopefully, physical insight, we often interchange the order of summations of 
infinite series (implemented as appropriate limits when necessary) with 
integration or differentiation, without asking whether the interchange is 
rigorously warranted from a mathematical viewpoint. The final results are 
then checked a posteriori for physical sense and consistency, and only at 
that stage, in doubtful cases, can we reassess our mathematics. In their own 
right, the ideas of the Dirac delta function, the Heaviside function, etc. — 
which have been widely used since the beginning of this century without a 
rigorous mathematical justification — reflect very well this ‘practical’ attitude. 

Before getting a bit more involved with the mathematical details, let us 
consider briefly the standard (practical) approach. 
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Intuitively, the Dirac delta function S(t) (we write it as a function of time 
only for future convenience) is defined to be zero when t^O and infinite at 
t = 0 in such a way that the area under it is unity. So, if a and b are two 
positive numbers, we can write 

S(t) = 0 t + 0 (2.59a) 

and 



" b 

mdt = 

J -a 



j 6(t}dt = 1 

J— OO 



(2.59b) 



which evidently implies 6(t)dt = \ b a 8{t)dt = 0. The weakness of the 
definition given by eqs (2.59a and b) is evident: strictly speaking, one can 
object that, first, 8(t) is not a function and, second, that eqs (2.59a) and (2.59b) 
are incompatible (if a function is zero everywhere except at one point, its 
integral, no matter what definition of integral is used, is necessarily zero). 

Nonetheless, one generally ignores these objections and proceeds by stating 
that the definition of eqs (2.59a and b) can be replaced by the following and 
equivalent definition: let g(t) be any well-behaved function which we may 
call a ‘test function’, then 

| g(tmdt = g{ 0) (2.60) 



the rationale behind eq (2.60) being that, by virtue of (2.59a), there is no 
contribution to the integral from anywhere where g(f) / g(0) so that g(t) can 
be replaced by the constant g( 0) which, in turn, can be moved outside the 
integral sign. The equality in (2.60) then follows from eq. (2.59b). 

Equation (2.60) is the usual definition commonly found in textbooks, 
which is sometimes called the ‘weak definition’ of 8(t) because it leaves 
undefined the value of the integral when one of the limits of integration is 
precisely zero. As a matter of fact, eq (2.60) and its consequences are 
compatible with many different assignments of the integrals 




(2.61) 



which are only subject to the condition J ^ 5(t)dt + f 8(t)dt — 1. 
Therefore, any definition that assigns values to the integrals of eqs (2.61 ) can 
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be called a ‘strong’ definition of 8(t). In practice, however, provided that no 
integration stops at zero (as is often the case), no strong definition is necessary 
and we can simply ignore the difference. On the other hand, when one of the 
integration limits is zero, we definitely need a strong definition, also in the light 
of the fact that there exist equations that are ‘weakly true but strongly false’. 

On physical grounds, the delta function may arise, for example, in a situation 
in which we consider a sudden impulsive blow applied to a mass m. If this 
impulsive force f(t) lasts from t=t 0 to t-t u from Newton’s second law we get 



r'l 


dv 


f{t)dt = 


m — ■ dt = 


b J 


1 dt . 



mdv = m{vi 



wo) 



( 2 . 62 ) 



stating that the impulse of f(t ) equals the change in momentum of the mass. 
Also, note that the precise shape of f(t) is irrelevant, the relevant quantity 
being its area between the two instants of time. So, let this area be unity and 
let, for simplicity, the initial velocity of the mass be vo = 0. If t = fo is very 
small we can simply ignore the motion of the mass during this time and 
assume that its momentum changed abruptly from zero to mv v As t\ — to 
gets smaller and smaller, the graph of momentum versus time approaches a 
step Heaviside function and — since (see eq (2.62)) the function f(t) is the 
slope of this graph — we are practically requiring f(t) to be the derivative of 
such a step function, i.e. infinite at t=t 0 and zero otherwise. No ordinary 
function possesses these properties; however, since we are not so much 
interested in f(t) itself but in the effect it produces, it can be convenient to 
introduce an entity — identified by the symbol 8(t — to) — that represents this 
infinitely high and infinitely narrow ‘unnatural’ force. In this light, 8(t) can 
physically represent a sudden blow applied at t o =0. 

So, returning to the general discussion, it is not difficult to show that the 
following (weakly true) properties of 8(t) hold: 



1. 



■b 



g[t)S(t - c)dt = 




a<c<b 

otherwise 



(2.63) 



Also 

„ 1 

2. 8{at) = S(t) 



(2.64) 



where a / 0. From the evenness of 8(t) (i.e. 5(— t) = <5{f) or 6(t — t) = S(t — f)) 
it follows that 



3. 




r)6{t — r)dr = 



g{T)S(T - t)dr = g{t) 

J-oc 



(2.65) 
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with the consequence that g(t ) * 8(t) = g{t ), meaning that 8(t) is the identity 
element in the operation of convolution. Next, if we consider the step 
Heaviside function 



m = 



0 *<0 
1 t> 0 



we can formally write 



4. 6(t) = -6(t) 
dt 



and also 

5. 0(t) = f S(r)dT 

J— 00 



( 2 . 66 ) 



(2.67a) 



(2.67b) 



Then, if we turn our attention to the derivative 8'(t) of 8(t), we may note 
that this function is not directly determined by the definition of 8(t) so that, 
consequently, it is meaningless to ask what this derivative is. Instead, in 
order for the ordinary rules of integration by parts remain valid, the approach 
is to define 8'(t) in such a way that 8’{t) = 0 when t^ 0, In symbols, if g(t) is 
a differentiable test function, we get 



r b d6(t) 

g(t) — 
- a at 



dt = g{t)8{t)\ b - 




dg(t) 



dt 



t = o 



and hence 

6. [ g[t)8’{t)dt= -g‘( 0) (2.68) 

so that, by extension, the nth derivative d ,n] (t) is required to vanish at 0, 
and for a well-behaved function g(t) we get 

7. f g{t)8 M {t)dt = ( - 1)"£ W (0) (2.69) 

J — Xr 

The formal ‘proofs’ of all the above properties are not difficult and are 
left to the reader. Other properties will be considered if and whenever needed 
in the course of future chapters. 



Copyright © 2003 Taylor & Francis Group LLC 




